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Abstract - The earlier-developed ab initio model and the kinetic Monte Carlo 
method (KMCM) are used to simulate precipitation in a number of iron-copper 
alloys with different copper concentrations x and temperatures T. The same simu- 
lations are also made using the improved version of the earlier-suggested stochastic 
statistical method (SSM). The results obtained enable us to make a number of 
general conclusions about the dependencies of the decomposition kinetics in Fe-Cu 
alloys on x and T. We also show that the SSM describes the precipitation kinet- 
ics in a fair agreement with the KMCM, and employing the SSM in conjunction 
with the KMCM enables us to extend the KMC simulations to the longer evolution 
times. The results of simulations seem to agree with available experimental data for 
Fe-Cu alloys within statistical errors of simulations and the scatter of experimental 
results. Comparison of results of simulations to experiments for some multicompo- 
nent Fe-Cu-based alloys enables us to make certain conclusions about the influence 
of alloying elements in these alloys on the precipitation kinetics at different stages 
of evolution. 

1 INTRODUCTION 

Studies of the composition and temperature dependencies of the precipitation kinetics 
in alloys, in particular, in multicomponent Fe-Cu-based steels used in many industrial 
applications, attract great attention PQ-[H]. For such studies, it seems to be useful to 
have a reliable information about the similar dependencies for more simple, binary alloys 
treated as the reference systems. For example, discussing the precipitation kinetics in 
two Fe-Cu-based steels, NUCu-140 and NUCu-170, Kolli et al. E] compared their 
experimental results to those for several binaries Fe-xCu with x between 1.34 and 1.5 
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at.%, while the two steels considered had notably different content of copper, 1.17 and 
1.82%, respectively. As many characteristics of nucleation and growth, in particular, sizes 
and density of precipitates, strongly vary with x = xq u , such comparison can be not quite 
adequate. The precipitation characteristics can also significantly depend on temperature, 
while measurements of such dependencies, particularly for the nucleation stage, often 
meet difficulties [5J E] . 

To obtain the quantitative information about the precipitation kinetics, particularly 
in the course of nucleation and growth, one can use simulations of these processes, if both 
the microscopic model and the methods of simulations used can be considered as reliable. 
For the Fe-Cu alloys, such reliable ab initio model has been developed by Soisson and Fu 
[H], and their detailed kinetic Monte Carlo simulations of nucleation and growth in the 
Fe-1.34Cu alloy at T=773 K revealed a good agreement with the available experimental 
data. 

Therefore, the first aim of this work is to use the model and the kinetic Monte Carlo 
method (KMCM) developed in [H] to study the concentration and temperature dependen- 
cies of the precipitation kinetics in binary Fe-Cu alloys for a number of copper concen- 
trations x and temperatures T, including those used by Kolli et al. [6] for the NUCu 
steels. Comparison of results of these simulations to the available experimental data en- 
ables us to assess reliability of the Soisson and Fu model at different x and T, and also to 
discuss the differences in the precipitation kinetics between the Fe-xCu binaries and the 
multicomponent Fe-Cu-based alloys with the same x = Xcu and T. 

The second aim of this work is to discuss possible applications of the earlier-suggested 
stochastic statistical method (SSM) [13] for simulations of precipitation kinetics in those 
cases when the KMCM meets difficulties. Such difficulties arise, for example, in treatments 
of the coarsening stage or at relatively low temperatures T < 300° C (typical for service of 
many nuclear reactors) when the KMC simulations become time-consuming [HldD]- Unlike 
the KMCM, the SSM allows parallelizing computer codes which can greatly accelerate 
computations, and this method seems also to be suitable for various generalizations, for 
example, for considerations of lattice misfit effects. However, in treatments of nucleation, 
some oversimplified models have been used in Ref. [13] which resulted in some fictitious 
breaks in simulated temporal dependencies. In this work we describe the improved version 
of the SSM free from these shortcomings and show that this version describes the main 
characteristics of nucleation, including the density and sizes of precipitates, in a good 
agreement with the KMCM. We also show that employing the SSM in conjunction with 
the KMCM enables us to extend the KMC simulations to the first stages of coarsening. 

In Sec. 2 we briefly discuss the model, the alloy states and the methods used in our 
simulations. In Sec. 3 we describe the improved version of the SSM employed in this work. 
The results of our simulations are discussed and compared to the available experimental 
data in Sec. 4. The main conclusions are summarized in Sec. 5. 



2 MODELS AND METHODS OF SIMULATIONS 



For our simulations we use the ab initio model of Fe-Cu alloys developed by Soisson and 
Fu and described in detail in Refs. [9j [TOj H3] . Here we only note that this model uses the 
following values of the binding energy between two copper atoms and between a copper 
atom and a vacancy, E^ Cu and E^ uv , for the n-th neighbors (in eV): 
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The high values of -EcuCu correspond to the strong thermodynamic driving force for pre- 
cipitation, while the strong attraction of a vacancy to copper atoms results in the strong 
vacancy trapping by copper precipitates discussed in detail in j9]. 

TABLE 1. The Fe-a;Cu alloy states considered and parameters of "thermodynamic" 
critical embryos for these states calculated by method of Ref. [T4] . 



Alloy state 


T K 


x, at% 


s 


Fr/T 


ATth 

1 v c 


R^, nm 


A 


773 


1.17 


0.247 


5.58 


15.4 


0.353 


B 


773 


1.34 


0.285 


4.38 


14.3 


0.344 


C 


773 


1.82 


0.393 


2.30 


12.2 


0.326 


D 


713 


1.34 


0.352 


2.47 


10.7 


0.312 


E 


663 


1.34 


0.425 


1.36 


8.2 


0.286 


F 


561 


0.78 


0.387 


1.11 


5.3 


0.247 


G 


873 


1.15 


0.163 


12.6 


28.3 


0.433 



The alloy states (below for short: states) used in our simulations are listed in Table 
1. The degree of supersaturation for each of these states is characterized by the reduced 
supersaturation parameter s introduced in Refs. [TBI [14]: 

s(x, T) = [x - x h (T)\/[x s (T) - x b (T)} (2) 

where the lower index s or 6 corresponds to the spinodal or the binodal. Values s < 1 
correspond to the nucleation and growth evolution type, and s > 1, to spinodal decom- 
position. 

The states A and C in Table 1 have the same temperature T and the copper content x 
as the above-mentioned steels NUCu-140 and NUCu-170 studied by Kolli et al. [21 E]. For 
the states B and G, kinetics of precipitation under thermal aging was investigated in Refs. 
[D El [31 [7] , while for the state F, precipitation under neutron irradiation was studied by 
Miller et al. [3]. For the states D and E, experimental data about precipitation kinetics 
are unknown to us, and our simulations are made to study the temperature dependence 
of this kinetics. 

In Table 1 we also present some parameters of the "thermodynamic" critical embryos 
for the states considered: the nucleation barrier F c , the total number of copper atoms 
within the embryo, N* h , and the critical radius R 1 ^ defined as the radius of the sphere 
having the same volume as Nl h copper atoms in the BCC lattice of a-iron with the lattice 
constant a=0.288 nm: 

Rf = a (3A"f787r) 1/3 = 0.142 (Af ) 1/3 nm. (3) 

The characteristics of critical embryos shown in Table 1 and Fig. 1 are calculated by 
the statistical method of Dobretsov and Vaks [J4] with the use of the ab initio model 
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Figure 1: (color online [II])- Concentration profiles Ac(R) = c(R)—x in thermodynamic 
critical embryos for the alloy states considered, where c(R) = cc u (r) is the mean copper 
concentration at the distance R from the center of embryo. 

by Soisson and Fu [5] mentioned above and the pair cluster approximation which is typ- 
ically highly accurate, particularly for dilute alloys, as both analytical studies [15] and 
comparisons with Monte Carlo simulations [16] show. Table 1 and Fig. 1 illustrate, 
in particular, the decreasing of the nucleation barrier F c and the embryo size Nf 1 with 
increasing supersaturation s. Table 1 also shows that for the given supersaturation s, 
lowering temperature T affects the critical embryo characteristics stronger than increas- 
ing concentration x [13J. As discussed below, the A*' 1 and Rf 1 values presented in Table 
1 are usually close to those estimated in KMC simulations and in experiments. 

Methods of both KMC and SSM simulations used in this work are described in detail 
in Refs. [9j [13]. Here we only mention some details of computations. For both the KMC 
and the SSM simulations we usually employed the simulation volume V s = (64a) 3 ; for 
the state G with a low supersaturation we used the bigger V s = (128a) 3 . In our KMC 
simulations we usually followed N^ MC ~ 10 13 Monte Carlo steps, which at V s = (64a) 3 
took about five months on a standard workstation. The SSM simulations for V s = (64a) 3 
took usually about a month on a standard workstation (with no parallelization of codes 
made) . 

3 IMPROVEMENTS OF STOCHASTIC STATISTI- 
CAL METHOD 

The original version of the SSM is described in Ref. [32]. To explain its improvements 
made in this work we first present the necessary relations from [13]. Evolution of a binary 
alloy is described by a stochastic kinetic equation written in the finite difference form for 
a short time interval St: 

5a = d(t + St) - a{t) = Scf+Y^ <H (4) 

where q is the occupation of site i by a minority (copper) atom averaged over some 
locally equilibrated vicinity of this site, and the "diffusional" term Scf corresponds to the 
average atomic transfer for a certain effective direct atomic exchange (DAE) model: 

Scf{c k } = M H 2 sinWA; - Ai)/2] St. (5) 

jnn ii) 

Here Aj is the local chemical potential of site % given by Eq. (20) in [13] . and M^- is the 
generalized mobility: 

^ = 7: / /^exp[/3(A 4 + A J )/2]. (6) 

The factor 6^ in (jBJ) is some function of local concentrations q given by Eq. (33) in [13] , 
while is an effective rate of exchanges between neighboring minority and host atoms, 
a and h (copper and iron), which is proportional to the product of analogous rates 7 QV 
and 7h v describing exchanges between a vacancy and an atom a and between a vacancy 
and an atom h, respectively: 

lf h = l a ,lh,v(t). (7) 
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Here the function u(i) defined by Eq. (36) in [13] determines the "rescaling of time" 
between the original vacancy-mediated exchange model and the effective DAE model 
used for simulations. Temporal evolution of this DAE model is described by the "reduced 
time" t r having the meaning of a mean number of atomic exchanges a ±=t h for the 
physical time t. The reduced time t r is related to this time t by the following differential 
or integral relations: 

dtr = <y%dt, 7^ = 7«vTfcvK*) J *= / K/lfhiK)- (8) 

J 

The form of the function t(t r ) is discussed below. As mentioned in Ref. [32], Eqs. 
dH)-(jHl) can be derived not only for a simplified "quasi-equilibrium" model used in [13] 
that disregards renormalizations of effective interactions in a nonequilibrium alloy (whose 
importance for the diffusion kinetics was noted by Nastar et al. [19]) but also with allowing 
for these renormalizations. 

The last term 5n{- in Eq. (pTJ is the fluctuative atomic transfer through the bond ij 
described by the Langevin-noise-type method: each 5n{ 3 is treated as a random quantity 
with the Gaussian probability distribution: 

W{8n{ 3 ) = A l3 exp[(-54) 2 /2Ai] (9) 

where Ay is the normalization constant. The dispersion Dij is related to the mobility 
Mij and the time interval 5t in Eq. (j5J) by the "fluctuation-dissipation" type relation: 

D t] = ((5n{ j ) 2 ) = 2M tJ 5t. (10) 

As discussed in detail in JT3] , for the non-equilibrium statistical systems under considera- 
tion, the Langevin-noise-type equations fj4l)- (fT0l) should be supplemented by the "filtration 
of noise" procedure eliminating the short-wave contributions to fluctuations Sn{ 3 as these 
contributions are already included in the diffusional term Scf obtained by the statistical 
averaging over these short-wave fluctuations. Therefore, in the last term of Eq. 01]), the 
full fluctuative transfer 8n io - should be replaced by its long-wave part Sn{j. It can be 
made by introducing a proper cut-off factor F c (k) in the Fourier- component 5nf a (k.) of 
the full fluctuation 8n i3 - = <5n£(R SQ ,) where H sa denotes the position of the ij bond 
center in the appropriate crystal sublattice a formed by these centers [T2] : 

Snt -»■ 5n fc (R sa ) = ^exp(-ikR sct ) 5n f (k) F c (k) 

k 

5n f (k) = 1 ^exp(zkR s ) 5n f (R sa ) (11) 

where N is the total number of lattice sites (or atoms) in the crystal. The cut-off factor 
F c (k) for the BCC lattice can be taken in the Gaussian-like form: 

F c BCC (k) = exp [-4# 2 (1 - cosy?icosv9 2 cosy?3)] (12) 

where (p v = k u a/2; k v is the vector k component along the main crystal axis v\ and a is 
the BCC lattice constant. At large g 2 3> 1, the expression (|T2|) is reduced to a Gaussian 
exp (— k 2 l 2 /2) with / = ga. Therefore, the reduced length g = l/a characterizes the 
mean size of locally equilibrated subsystems. 
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Figure 2: (color online [17J). Evolution of the free energy per copper atom, F(t r ) (a), 
and the density d p (t r ) of supercritical precipitates or their total number N p (t r ) within 
simulation volume V s = (64a) 3 (6), obtained for the state A in the SSM simulations with 
different g and t . Curves 1, 2, and 3 (red, green and blue online [T7]) correspond to 
t = 40 and g =1.4, 1.35, and 1.45, respectively. Curves 4 and 5 (red on-line) correspond 
to g Q = 1.4 and t =50 and 30, respectively. Thin vertical line corresponds to t rN = 18. 

This size, generally, varies with the aging time t or t r . In particular, after creation of 
a supercritical precipitate, the degree of local equilibrium in the adjacent region should 
increase with respect to other regions where such precipitates are not born yet. Therefore, 
after completion of nucleation at some reduced time t r ^ (which can be estimated as the 
time for which the density of precipitates reaches its maximum), the alloy should rapidly 
approach the two-phase equilibrium, and the length I = ga should become large, which 
should lead to a sharp decrease of the ffuctuative terms 8n* = oW c in Eqs. ( TTTj) and 

To describe this physical picture with the minimal number of model parameters, the 
time dependence g(t r ) was approximated in [13] by a simple one-parametric expression 
(71) with a break at t r = t r M- In the present work we use for g(t r ) a smooth two- 
parametric expression: 

g(t r )=g (l + t 2 r /t 2 ) (13) 

where go and t Q have the order of magnitude of the reduced critical radius R c /a and the 
reduced nucleation time t r ^, respectively. 

TABLE 2. Parameters go and to in (TT3]) . reduced nucleation time t r N, and maximum 
precipitate density d max (in 10 24 m~ 3 ) for the alloy states considered. 



Alloy state 


% 


to 


Un 


jSSM 
^max 


jKMC 
max 


A 


1.4 


40 


~ 20 


3.4 


4.0 


B 


1.4 


30 


~ 20 


5.4, 5.8 


6.1 


C 


1.45 


20 


~ 10 


14.9 


14.2 


D 


1.4 


10 


~ 10 


9.7 


10.1 


E 


1.6 


10 


~ 10 


13.6 


11.2 


F 


1.6 


5 


~ 5 


9.9 


8 



To estimate values go and t within the SSM, we can use the "maximum thermo- 
dynamic gain" principle discussed in detail in [13]: we suppose that the "most realistic" 
values of these parameters correspond to the minimum of the free energy of an alloy after 
completion of nucleation. To illustrate this variational method of estimating go and to, 
in Fig. 2 we present the SSM-simulated temporal dependencies of the free energy per 
copper atom for the state A in Table 1 at several go and to close to their "optimal" 
values. The free energy F = F(ci) was calculated according to Eq. (24) in [13] with the 
simulated Cj = Ci(t r ) dependencies. For simplicity, the initial state was taken uniform: 
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Figure 3: (color online [E]). Evolution of the density of precipitates, d p (t r ) (a), and their 
mean radius R(t r ) (b) obtained in the SSM simulations for the state A in Table 1. Dotted 
vertical lines correspond to the values ti r in Eqs. ffT6j) . 

Cj(0) = c=const, thus the initial increase of F at t r < 0.3 t r ^ seen in Fig. 2a is related 
just to switching-on fluctuations at t r = 0. Fig. 2a shows that the influence of this 
spurious increase of F ceases only at t r > 2t r ^, while at too long t r > 4t r N, the fluc- 
tuations are effectively switched-off according to Eqs. (TTTjl -l TTB"]) . Therefore, to estimate 
parameters go and to, we usually consider the interval 2t r N ^t r > 4t r N. For the most of 
states considered, the free energy F(t r ) has here a distinct minimum with respect to g 
and t , as illustrated by Fig. 2a. Fig. 2b illustrates sensitivity of the precipitate density 
dp to the choice of go and to- The resulting scatter in simulated d p (t r ) usually lies within 
statistical errors of simulations. 

In Table 2 we present the values of go and to in ([13]) estimated as described above. 
For the state G the simulations are time-consuming, thus such estimates have not been 
made in this work. To assess adequacy of our SSM simulations, in the two last columns of 
Table 2 we compare the values of the maximum density d max of supercritical precipitates 
obtained in the SSM simulations to those obtained in the KMC simulations. For the state 
B in Table 2, two values d^Jf correspond to the two simulations (shown in Fig. 6) with 
the different random number sets. We see that the SSM and the KMC results for d max 
agree within statistical errors of simulations. 

The precipitation kinetics is usually characterized by the density and the mean radius 
of supercritical precipitates, d p {t) and R(t), defined by the following relations: 

d P (t) =J2"k(t)/V s , 

k 

R(t) = J2^(t)R k /Y,^(t)- (14) 

k k 

Here Vk(t) is the number of clusters containing k copper atoms, V s is the simulation vol- 
ume, Rk is defined similarly to Eq. ([3]): Rk = a(3/c/87r) 1 / 3 , and sums over k in (fbil) include 
only clusters with k > N c where N c is the "critical" size chosen. As discussed in [HI H3] 
and below, the exact choice of this size (if reasonable) usually does not significantly affect 
the d p (t) and R(t) values in (TH1) . Therefore, we take N c close to its "thermodynamic" 
value Nf in Table 1 putting JV C =15, 15, 12, 11, 8, 5 and 28 copper atoms for the state 
A, B, C, D, E, F and G, respectively. Temporal dependencies d p (t r ) and R(t r ) obtained 
in our SSM simulations with such N c are illustrated by Fig. 3 for the state A; for the 
rest states considered, these dependencies are similar. Fig. 3 illustrates the sequence 
of four well-defined stages of precipitation accepted in the classical theory of nucleation 
(see, e. g., [20]): incubation, nucleation, growth and coarsening, as well as the presence of 
an "intermediate" stage between growth and coarsening (seen also in the simulations of 
Soisson and Martin |20j and discussed in detail by Barashev et al. [21]) which corresponds 
to the beginning of Ostwald ripening when the bigger precipitates already start to grow 
due to dissolving the smaller ones but the latter do not disappear yet. For brevity, this 
intermediate stage will be called "pre-coarsening". 

Let us now discuss the "rescaling of time" function t(t r ) in Eq. ([8]) determined by the 
temporal dependence of the effective direct exchange rate 7^ = 7cuF e (^)- F° r simplicity, 
this dependence was approximated in [T3] by a simple two-parametric expression (77) 
which included sharp breaks at some t r . More adequate interpolations for t{t r ) can be 
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obtained from comparison of the SSM and the KMC results for the density and the mean 
size of precipitates, that is, from an approximate solution of two equations: 

dp MC (t) = dp SM (t r ), RKMo{t) = RssM(t r )- (15) 

Our estimates of functions t(t r ) for which both equations ( 115]) are satisfied with a reason- 
able accuracy showed that for all alloy states considered, these functions have a similar 
form characterized by the presence of four intervals of t r with an approximately constant 
value of the derivative D = dhxt/dhxt r oc l/7cuF e within each interval. These intervals 
(i) are determined by the inequalities: 

(1) t r < t r i, (2) t r i < t r < t r 2, 

(3) t r2 < t r < tr3, (4) t r3 < U, (16) 

where the point t r i, t r 2 or t r3 approximately separates the stages of nucleation and growth, 
growth and pre-coarsening, or pre-coarsening and coarsening, respectively, as illustrated 
by Fig. 3. Therefore, within each interval (i) we use a simple power interpolation 

t = Ait®* (17) 

with the values of parameters t ri , Di and A4 given in Table 3. For the state F with a 
relatively low T = 561 K, our KMC simulations (illustrated by Fig. 11 below) are time 
consuming and reach only the growth stage, thus the interpolation (IT7|) has not been 
constructed for this state. 

Functions t(t r ) obtained are shown in Fig. 4. To compare the precipitation parameters 
at different temperatures T for which the equilibrium vacancy concentration c^ 9 (T) can 
be very different, in Fig. 4 we use the "scaled" time t s (t r ) defined as 

t s (t r ) = t(t r )c e J(T)/C(773K), (18) 

to cancel the standard scaling factor 1/c^ 9 in the t s value. For the c e ^(T) we use the 
Soisson and Fu estimate [9]: 

C(T) = exp (-4°7n e{ OT = 2.18 eV. (19) 

The SSM values dp SM (t) and RssM(t) found using Eqs. (TT7T) and Table 3 are compared to 
the appropriate KMC results in Figs. 5-10. Note that the interpolation ([17]) includes some 
minor breaks at points t r = t r i which are seen in Fig. 4. These breaks can be removed 
by introducing some smooth matching of two lines (1171) adjacent to each point t ri , which 
also leads to the better agreement between the SSM and the KMC results in Figs. 5-10. 
However, to avoid introducing too many interpolation parameters, in this work we use a 
more simple interpolation (TIT]) which seems to provide a sufficiently adequate description 
of precipitation. 

As discussed in Ref. [13], the significant changes of derivatives Di ~ V(7cuFe)i 
tween different intervals i seen in Fig. 4 can be related to the strong vacancy trapping by 
copper precipitates and to the changes of scale of this trapping in the course of precipita- 
tion. In particular, the relatively low values of derivatives D% for the pre-coarsening stage 
can be related to the strong vacancy trapping for this stage (illustrated by Fig. 6a in [9]) 
which should result in a sharp increase of the effective rate 7cup e i n Eq. (jHJ) P3]. Note 
also similarity of the curves t s (t r ) for different alloy states in Fig. 4, as well as rather 
smooth variations of parameters Di and Aj in Table 3 under changes of concentration 
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Figure 4: (color online [17]). The scaled physical time t s (t r ) defined by Eqs. (TTTj) and 
ffTB"]) versus the reduced time t r . Each curve corresponds to the alloy state indicated by a 
symbol near this curve. 

x and temperature T corresponding to these different states. It seems to reflect a great 
similarity of the vacancy trapping effects within each interval i considered for the different 
alloy states. This similarity can be used for various SSM-based extrapolations of KMC 
simulations, in particular, for the SSM-based simulations of precipitation in Fe-Cu alloys 
at most different concentrations x and temperatures T with the use for the parameters 
t ri , Di and in Eqs. ( I16p and ( I17p some interpolations between their values presented 
in Table 3. This similarity is also used for the extrapolations of KMC simulations to the 
first stages of coarsening discussed below. 



TABLE 3. Values of parameters t r i, Di and Ai (in hours) in Eqs. ( TT6|) and (I17p . 



Alloy state 




t T 2 


tr3 


D x 


D 2 


£3 


D 4 


A! 


A 2 


^3 


A, 


A 


18 


506 


5500 


0.7 


0.23 


0.04 


0.3 


0.061 


0.24 


0.77 


0.08 


B 


7.4 


192 


7240 


1.1 


0.33 


0.11 


0.4 


0.021 


0.097 


0.30 


0.023 


C 


2.9 


70 


3000 


1.85 


0.23 


0.12 


0.35 


0.013 


0.071 


0.11 


0.018 


D 


5 


60 


2300 


1 


0.23 


0.08 


0.27 


0.027 


0.14 


0.26 


0.026 


E 


3.2 


60 


2200 


1.5 


0.19 


0.08 


0.3 


0.015 


0.16 


0.25 


0.012 



To conclude this section we note that the simplified treatment of fluctuations based 
on Eqs. (l9|)- (fT3|) can adequately describe nucleation and growth only when these two 
processes are sufficiently separated from each other, so that switching off fluctuations 
after completion of nucleation implied by Eq. ( Fl3|) can be appropriate. Figs. 5-10 show 
that this condition is more or less satisfied for the first five states in Table 1. At the 
same time, for the states F and G, that is, at low temperatures T < 300° C or at low 
supersaturations s < 0.2, the nucleation, growth and pre-coarsening stages overlap very 
strongly, as Figs. 11, 12 and 17 show, and the simple model fl9|)-f JT3|) implying the type of 
evolution shown in Fig. 3 can hardly be adequate. Therefore, for the states F and G, only 
KMC simulations are presented in the present work. Further refinements of the SSM are 
evidently needed to employ this method at the low temperatures or low supersaturations 
mentioned. 

4 RESULTS OF SIMULATIONS AND COMPARI- 
SON WITH EXPERIMENTS 

The results of our simulations together with the available experimental data for some 
Fe-Cu and Fe-Cu-based alloys are presented in Figs. 5-17. Figs. 5-12 show the temporal 
evolution of the density and the mean radius of precipitates, d p (t) and R(t). Note that for 
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Figure 5: (color online [17J). Evolution of the density of precipitates, d p (t) (top figure), 
and their mean radius R(t) (bottom figure) obtained in our simulations for the state A with 
2=1.17 (in at.%, here and below), T=773 K for the critical size N c = 15 copper atoms. 
Right scale shows the number of precipitates within simulation volume V s = (64a) 3 . 
Solid curves correspond to the KMCM, and dashed curves, to the SSM. Points show 
experimental data by Kolli et al. [6] for the multicomponent steel NUCu-140 with the 
same xcu=l-17 and T=773 K for N c = 11. 

Figure 6: (color online [17]). The same as in Fig. 5 but for the state B with x=1.34, 
T = 773 K for N c = 15. Two dashed curves (red and green on-line) correspond to the two 
SSM simulations with different random number sets. Symbols correspond to experimental 
data for the state B: circles, to [1J; triangles, to [2]; and rhombus, to [3J. 

the state F in Fig. 11, this evolution is described in terms of the scaled time t s defined by 
Eq. (JTSjl (with replacing t r —>•£), as the equilibrium vacancy concentration c e v q (T) at the 
low T=561 K considered is negligibly small, while actually (in particular, in experiments 
[3]) the precipitation at such low T occurs only in irradiated materials where the vacancies 
(together with the interstitial atoms) are formed due to irradiation. 

Figs. 13-16 illustrate the concentration and temperature dependencies of the max- 
imum precipitate density d™ ax = d max , as well as of some temporal characteristics of 
precipitation, t max , t c and t c,0A , defined by the relations: 

d p (t max ) = d max , (20) 
t c = t(t r3 ), (21) 
d p (t c >°- 1 )=0.ld max (22) 

where the reduced time t r3 and the function t(t r ) in ( 121 j) are the same as in Eqs. (TT6"j) 
and (I17p. The time t max usually corresponds to the completion of nucleation or beginning 
of growth; the time t c approximately corresponds to the onset of coarsening, and t c ' 01 
can characterize the time of completion of the first stage of coarsening and beginning 
of its more advanced stages. For brevity, t max will be called the "nucleation time", and 
t c ' 01 , the "advanced coarsening time". In Figs. 13-16 we also present the experimental 
estimates of t max and t c ' 01 but not t c as such estimates are usually not certain for the 
onset of coarsening. 

Fig. 17 illustrates temporal evolution of the precipitate size distributions observed in 
our simulations. In caption to this figure we use the times t ,a and t c,a defined similarly 
to the "nucleation time" t max and the "advanced coarsening time" t c,ai in Eqs. ( 12 Op and 

d p (t N > a ) = ad max , t N > a < t max , (23) 
d p (t c ' a ) = ad maxi t c < a > t max (24) 

where the number a is less than unity. Times t N,a qualitatively correspond to the nucle- 
ation stage, and t c,a , to the coarsening stage. 

Let us discuss the results presented in Figs. 5-17. First, Figs. 5-12 and 17 show that 
the above-mentioned strong attraction between a vacancy and a copper atom leads to a 
great difference in the precipitation kinetics between iron-copper alloys and alloys with no 
such attraction, such as the alloys described by simplified models with time-independent 
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Figure 7: (color online [17J). The same as in Fig. 6 but for seven different KMC simula- 
tions for N c = 10. Thick solid line (red on-line) corresponds to the same KMC simulation 
as that shown in Fig. 6 by thick line for the choice N c = 15. 

Figure 8: (color online [17]). The same as in Fig. 5 but for the state C with £=1.82, 
T = 773 K for iV c = 12. Points show experimental data by Kolli et al. [5] E] for the 
multicomponent steel NUCu-170 with the same xcu = 1-82, T = 773 K for N c = 11. 

effective direct exchange rates 7^ in (jSJ) for which this kinetics is illustrated by Fig. 3, 
or the models with a relatively weak vacancy-minority atom interaction for which this 
kinetics is illustrated by Fig. 1 in [20] or Fig. 4 in [21]. For such simplified models, 
the presence of five well-defined stages of evolution shown in Fig. 3 is characteristic, 
including the distinct pre-coarsening stage mentioned above. On the contrary, in Figs. 
5-12 such an intermediate pre-coarsening stage is not seen, the "pure nucleation" and 
"pure growth" stages are relatively short (if exist at all), and the nucleation, growth and 
pre-coarsening stages significantly overlap each other. It is also illustrated by the size 
distribution functions Vk{t) in Fig. 17 which seem to imply that the nucleation, growth 
and pre-coarsening processes at 0.5t max < t < t max occur simultaneously. As discussed 
in [13], these kinetic features seem to be mainly related to the strong vacancy trapping 
by copper precipitates, which leads to a great acceleration of both the growth and the 
Ostwald ripening processes as compared to the alloys with no such trapping. 

Figs. 5-10 also show that the simplified direct-atomic-exchange model fl4])-f lT2"]) us- 
ing the above-mentioned maximum thermodynamic gain principle and the simple inter- 
polation (|17p describes the precipitation kinetics in a fair agreement with the KMCM. 
Discrepancies between the SSM and the KMCM results lie usually within statistical er- 
rors of simulations, and these discrepancies can be still more reduced if the smoothed 
interpolations mentioned above are used instead of more simple expressions (fl7|) . 

Figs. 5, 6 and 8-10 also illustrate opportunities to use the SSM to extrapolate the 
KMC simulations to the longer evolution times. As mentioned, the KMC simulations of 
coarsening are time-consuming, unlike the SSM simulations. At the same time, the above- 
described considerations about the physical nature of the "rescaling of time" function 
t(t r ) in ffT7"|) enable us to expect that this function preserves its form for the first stages 
of coarsening at least up to t ~ t c,ai , until the later, more advanced stages of coarsening 
(including the asymptotic Lifshits-Slyozov- Wagner stage [221 123]) start. Therefore, we can 
use the SSM simulated d p (t r ), R(t r ) and Eq. f[T7|) to estimate the d p (t) and R(t) values 
for those t for which the KMC results are not available. This is shown in Figs. 5, 6 and 
8-10. 

Figs. 6 and 7 also illustrate some methodical points. Two dashed curves in Fig. 6 (red 
and green on-line) correspond to two SSM simulations with the different random number 
sets, thus their difference illustrates the statistical scatter of the SSM simulation results. 
We see that this scatter is significant only at the very end of simulations when the total 
precipitate number N p becomes small. Similarly, seven different KMC simulations shown 
in Fig. 7 illustrate the statistical scatter of the KMC simulation results. Influence of the 
choice of the "critical" size N c in Eqs. ([14]) on the simulated d p (t) and R(t) is illustrated 
by Figs. 6 and 7 which correspond to the same alloy state B but to the different iV c , 15 
and 11 copper atoms, respectively, while the thick solid curve in both Fig. 6 and Fig. 
7 (red on-line in Fig. 7) corresponds to the same KMC simulation. We see that the 
variations of N c used make noticeable effects on the simulated d p (t) and R(t) only at the 
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Figure 9: (color online [17J). The same as in Fig. 5 but for the state D with £=1.34, 
T = 713 K for iV c = 11. 

Figure 10: (color online [UJ). The same as in Fig. 5 but for the state E with x=1.34, 
T = 663 K for iV c = 8. 

first stages of nucleation, while later on such effects become insignificant. 

Let us discuss the concentration and temperature dependencies of the kinetic charac- 
teristics of precipitation presented in Figs. 13-16. For the maximum precipitate density 
d max , these dependencies shown in Figs. 13 and 15 are mainly determined by the reduced 
supersaturation s (defined in Eq. §Z§ and Table 1) which characterizes the scale of the 
thermodynamic driving force for precipitation. For the given temperature T or the given 
concentration x, the d max value increases with s, and at low s < 0.3 or high T > 800 K 
this rise is rather sharp. At the same time, Fig. 15 shows that at not low s > 0.3 — 0.4 and 
not high T, the d max value changes with T more slowly, and at T < 700 K, the further 
lowering temperature makes little effect on the d max . It may imply that the slowing down 
of kinetics due to the strong vacancy-copper atom correlations discussed below becomes 
important for these T. The temporal characteristics of precipitation shown in Figs. 14 
and 166 decrease with increasing supersaturation s, which can be explained by an increase 
of the thermodynamic driving force. At the same time, Fig. 16a shows that at not low 
s > 0.3 — 0.4 and not high T < 750 K, the scaled times t™ ax , t c s and t c s ' w vary with 
temperature rather weakly. It can be explained by an interplay between an increase with 
lowering T of both the thermodynamic driving forces which promote the evolution, and 
the vacancy-copper atom correlations which reduce the copper diffusivity Dc u [9] and 
thereby slow down the evolution. 

Let us discuss the precipitate size distributions v k (t) shown in Fig. 17. First, we 
see that these distributions are usually rather broad, and sizes of different precipitates 
are typically very different. Therefore, characterization of these sizes by only their mean 
value R(t) used in Figs. 5-12 is oversimplified and incomplete. It is true not only for 
the coarsening stage (for which a great difference in the precipitate sizes is natural as 
the bigger precipitates coarsen due to dissolving the smaller ones) but also for all other 
stages of precipitation. It can be related to the strong overlapping of the nucleation, 
growth and pre-coarsening stages mentioned above. Comparison of the size distributions 
u k (t) at t=t max and t=\.\t max (that is, frames G2 and G3, B2 and B3, or D2 and D3 in 
Fig. 17) also illustrates the significant overlapping of growth and early coarsening stages. 
We see that the great majority of precipitates at t=l.lt max still continue to grow due to 
the absorption of copper atoms from the matrix, even though some smallest precipitates 
already start to dissolve. Such significant overlapping of growth and early coarsening 
stages agrees with the observations by Mathon et al. for the Fe-1.34Cu alloys (the state 
B in Table 1) [3], and it was also noted by Kolli et al. for the NUCu steels [6]. Second, 
the comparison of the size distributions Vk{t) at t=t max and t=l.lt max also illustrates 
very sharp variations of these distributions with the evolution time t, in particular, at t 
corresponding to the beginning of coarsening, which is not clearly seen in the "averaged" 
description of Figs. 5-12. Third, Fig. 17 shows that the precipitate size distribution u k (t) 
for the state G with the highest temperature T=873 K is, generally, much more uniform 
than Vk{t) for the states B and D with the lower T=773 K and 713 K, particularly for 
the early coarsening stage, which is illustrated by frames G4, B4 and D4. In particular, 
in frames B4 and D4 we observe only one very big cluster with A;~600 copper atoms, 
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Figure 11: The same as in Fig. 5 but for the state F x=0.78, T = 561 K for N c = 5. 



Figure 12: The same as in Fig. 5 but for the state G with £=1.15, T = 873 K for N c = 28. 
Solid symbols show experimental data by Shabadi et al. [7J for the following alloys aged 
at T = 873 K: circles, binary Fe-1.15Cu; squares, ternary Fe-l.14Cu-0.99Mn. Right scale 
shows the number of precipitates, N p , within simulation volume V s = (128a) 3 . 

numerous small clusters with /c<200 copper atoms, and 3-4 "middle-sized" clusters with 
&;~200-300 atoms (the same features are also observed in Uk(t) for the state A not shown 
in Fig. 17 on considerations of space), in the great difference with frame G4 (as well as 
G2 and G3) where the cluster size distribution is rather uniform. This difference can be 
related to the weakening of the vacancy-copper atom correlations at high T which can 
enhance the copper diffusivity Dc n and thereby promote growth of many big precipitates 
for the state G unlike the states A, B and D, but these points need further studies. 
Finally, the lower row of Fig. 17 illustrates features of precipitation at low temperatures 
T ~ 300° C when the vacancies (necessary for the atomic diffusion) are provided by 
irradiation. Figs. F1-F3 illustrate a very strong overlapping of the nucleation growth and 
pre-coarsening stages, while the coarsening stage was not reached in these our simulations. 

Let us now compare the simulation results to the available experimental data. For 
the state B and the KMC simulations shown in Figs. 6 and 7 by thick lines, the detailed 
comparison with various experiments was given by Soisson and Fu [9] who concluded that 
the predictions of simulations are reliable. The results presented in Figs. 6, 7, 14 and 
16 can complement their discussion by the two points. First, Figs. 6 and 7 confirm that 
the disagreements between simulations and experimental observations seem usually to lie 
within statistical errors of simulations and the scatter of experimental results. Second, 
Figs. 6, 14 and 16 show that the SSM-based extrapolations of KMC simulations for the 
first stages of coarsening seem to agree with the observations as well. 

For the state G with x — 1.15, T — 873 K, our simulations are compared to the 
data by Shabadi et al. [7] in Fig. 12. As supersaturation s and the precipitate density 
dp for this state are rather low, these simulations are time-consuming and include only 
nucleation and growth stages, while the data by Shabadi et al. seem to correspond to the 
longer aging times and have significant errors. Within these errors, the simulation and 
experimental results in Fig. 12 can be considered as agreeing with each other, particularly 
for the maximum precipitate density d max , though the simulated evolution times can be 
somewhat shorter than the observed ones. Solid squares in Fig. 12 correspond to a 
ternary Fe-Cu-Mn alloy and illustrate the effect of the third alloying element Mn on the 
precipitation kinetics. The presence of Mn seems to lead to an increase of the maximum 
precipitate density d max by about twice with respect to the analogous Fe-Cu binary, in a 
qualitative contrast with the effect of alloying elements on the d max value in the NUCu 
steels discussed below. 

Figure 13: Concentration dependence of the maximum density of precipitates, d max , for 
decomposition of Fe-xCu and Fe-Cu-based alloys at T = 773 K. Open symbols here and 
below correspond to our simulations, and solid symbols, to experiments. Solid circle 
corresponds to experiments [2] for the state B, while solid square and solid triangle, to 
experiments [5, 6J for the multicomponent steel NUCu-140 with xc u =l-17 and NUCu-170 
with xcu=l-82, respectively. Dashed lines here and below are given to guide the eye. 
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Figure 14: Concentration dependence of temporal characteristics of precipitation: t max 
(circles), t c (triangles), and t c,ai (squares), defined in the text by Eqs. (l20i) . (l21i) and 
(12"2"|) . respectively, at T = 773 K. Solid circle and solid square at xc u =l-17 or x Cu =1.82 
correspond to the data ^ |BJ for the steel NUCu-140 or NUCu-170, respectively, while 
solid circle and solid square at £=1.34 correspond to the data [1]- [3] for the state B. 

Figure 15: Temperature dependence of the maximum density of precipitates in the course 
of decomposition of Fe-xCu and Fe-Cu-based alloys. Open circles correspond to the states 
B, D and E with £=1.34; open square, to the state A with £=1.17; and open triangle, to 
the state G with x=1.15. Solid circle corresponds to experiments [2] for the state B with 
x=1.34; solid square, to experiments [6J for the steel NUCu-140 with xcu=l-17; and solid 
triangle, to experiments [7j for the state G with £=1.15. 

In Figs. 5, 8, and 13-16 we compare the simulation results for the state A or C with 
T = 773 K and x = 1.17 or 1.82 to the data by Kolli et al. E] for the multicomponent 
steels NUCu-140 and NUCu-170 with the same T and icu- In addition to copper, these 
steels contain a number of alloying elements: C, Al, Ni, Si, Mn, Nb, P and S, 5.49% on 
the total in NUCu-140, and 5.83%, in NUCu-170, while the partial concentrations of each 
alloying element in these two steels are very close to each other [5j[6]. Therefore, differences 
in the precipitation kinetics for these two steels can be mainly related to the difference 
in the copper content xqu- Then comparison of this kinetics for each of these steels to 
that for the analogous binary Fe-xCu alloy can elucidate the effect of alloying elements on 
the precipitation at different xc u - Qualitatively, these problems were discussed by Kolli 
et al. [6J. Our simulations enable us to consider these points quantitatively. Let us also 
note that the critical sizes N c for these two steels estimated by Kolli et al. [SJ |BJ: N c ~11 
copper atoms, are rather close to the estimates for their binary analogues presented in 
Table 1: N?(A) ~ 15, Nf{C) ~ 12 copper atoms. 

Let us first discuss the nucleation and growth stages illustrated by Figs. 5 and 8. For 
the NUCu-140 or NUCu-170 steel this corresponds to t < 1 h and t < 0.25 h, respectively, 
and the maximum precipitate density d max in each steel is lower than that in its binary 
analogue by about three times. However, for the NUCu-140, both the values and the 
temporal dependencies of d p {t) and R(t) in Fig. 5 seem to not greatly differ from those 
simulated for the Fe-1.17Cu alloy, particularly for the nucleation stage, and the nucleation 
time t max can be similar, too. On the contrary, for the NUCu-170, the data at t=0.25 h 
shown in Fig. 8, particularly for the R(t) value, sharply disagree with those simulated for 
the Fe-1.82Cu alloy, while the nucleation time t max exceeds that for the Fe-1.82Cu alloy 
by an order of magnitude. 

Therefore, our comparison seems to imply that the effect of almost the same content of 
alloying elements on the nucleation kinetics in NUCu-170 with the higher copper content 
xcu=l-82 is much stronger and qualitatively different from that in NUCu-140 with the 
lower £cu=l-17. Physically, such conclusion does not seem to be natural. In this connec- 
tion we note that this conclusion is mainly based on the data for the mean precipitate 
size in NUCu-170 at t=0.25 h reported by Kolli and Seidman [S]: R — 1.2 nm. This value 
much exceeds the critical radius R c ~ 0.3 nm estimated for this steel in [5]: R ~ 4R C . It 
should imply that in the course of the nucleation stage (supposed in [5] for the NUCu-170 
at i=0.25 h to explain a steep rise of the precipitate density d p (t) between t=0.25 h and 
t=l h seen in Fig. 8) the new-born precipitates grow extremely fast. Such a very sharp 
growth at the early nucleation stage seems to be very unusual and, to our knowledge, 
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Figure 16: Temperature dependence of the same temporal characteristics of precipitation 
as those in Fig. 14, but in terms of the scaled time t s defined by Eq. (fl8l) : t™ ax (circles), 
t c s (triangles), and t^ ,ai (squares). In figure a, the states B, D and E have the same 
x=1.34, while solid circle and solid square correspond to experiments [2] for the state B 
with 2=1.34. In figure b, the state A corresponds to £=1.17; solid circle and solid square 
at T=773 K correspond to the data [6J for the steel NUCu-140 with xcu=l-17; and solid 
circle at T=873 K corresponds to experiments [7j for the state G with x=1.15. 

Figure 17: Numbers of clusters containing k copper atoms, observed in our KMC 

simulations. The first, second, third and fourth row corresponds to the state G, B, D 
and F in Table 1, respectively. Frame Gl, G2, G3, or G4 corresponds to the time t 
equal to t^' - 47 , t max , l.li mm ^tf s > '* r , or t c ' - 67 ; frame Bl, B2, B3, or B4, to t^' - 61 , t max , 
l.lf™*=t c > - 94 , or t c > om ; and frame Dl, D2, D3, or D4, to t N ' - 6 , t max , 1.1 t max =t c ' om , 
or £ c,a61 , respectively, where the time t N,a , t max or t c,a is defined by Eq. (|23|) . ( 120]) . or 
(I2~4j) . Frame Fl, F2, or F3 corresponds to the scaled time t s equal to tf' 0A8 , t^' - 94 , or 
t™ ax , where the time tf' a is defined by Eq. (1251) with replacing t— >-i s , while t™ ax = 1.45 h 
corresponds to the end of our simulations for the state F. 

was never observed in either experiments or simulations, as illustrated by Figs. 5-12. 
Therefore, the data about R(t) in NUCu-170 at t=0.25 h reported in [5] should possibly 
be taken with some caution. 

For the coarsening stage, the results presented in Figs. 5, 8 and 14 fully agree with 
the main conclusions of Kolli et al. [6] about a very strong slowing down of coarsening in 
the NUCu steels with respect to binaries Fe-Cu. In particular, the advanced coarsening 
time t c,0A for each of these steels exceeds that for its binary analogue by about two orders 
of magnitude. At the same time, the dependencies of this advanced coarsening time t c,0A 
on the copper content xqu for the NUCu steels and for their binary analogues shown in 
Fig. 14 seem to be similar. 

The strong slowing down of coarsening in a ternary Fe-Cu-Mn alloy with respect to its 
binary analogue Fe-Cu was also observed by Miller et al. under neutron irradiation [4|. 
On the contrary, the effects of alloying elements on the nucleation and growth kinetics in 
Fe-Cu-Mn alloys and in NUCu steels seem to differ qualitatively: according to Figs. 12 
and 13, the d max value for an Fe-Cu-Mn alloy is by about twice higher, while for each of 
NUCu steels, it is by about three times lower than that in its binary analogue. Therefore, 
for the coarsening stage, the effects of alloying elements on the decomposition kinetics in 
the multicomponent Fe-Cu-based alloys seem to be much more universal than those for 
the earlier stages of precipitation. 

Such a universal slowing down of coarsening in the multicomponent Fe-Cu-based alloys 
with respect to their binary analogues Fe-Cu can be related to a significant segregation of 
alloying elements on the surface of precipitates [l]-[n] which can reduce the surface energy 
and thereby the thermodynamic driving force for coarsening. It can also be related to 
a weakening of the vacancy trapping at surfaces of precipitates due to this segregation. 
However, quantitative estimates of these effects seem to be absent yet. 
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5 Conclusions 



Let us summarize the main results of this work. The earlier-developed ab initio model 
and both the kinetic Monte Carlo (KMC) and the stochastic statistical methods are 
used to simulate the precipitation kinetics for seven binary Fe-Cu alloys with different 
copper concentrations x and temperatures T. Comparison of results obtained to available 
experimental data and to other simulations enable us to make a number of conclusions 
about kinetic features of precipitation in both the binary Fe-Cu and the multicomponent 
Fe-Cu-based alloys. 

First, we find that due to the strong vacancy trapping by copper precipitates, the 
precipitation kinetics in iron-copper alloys for all x and T considered differs notably from 
that observed for the alloys with no such trapping: the "pure nucleation" and "pure 
growth" stages are relatively short, the nucleation, growth and coarsening stages signifi- 
cantly overlap, while the intermediate "pre-coarsening" stage observed in some simulations 
for simplified alloy models (illustrated by Fig. 3 of this work and by Fig. 1 in Ref. |2Uj ) 
is absent. In this connection, the presence of this pre-coarsening stage in simulations of 
precipitation in irradiated Fe-Cu alloys made by Barashev et al. [21] can be related just 
to some oversimplifications of their model. 

The concentration and temperature dependencies of the maximum precipitate density 
d max , the nucleation time t max , and the advanced coarsening time t c ' 01 defined by Eqs. 
(120]) and (|2"2"j) are illustrated by Figs. 13-16. At low supersaturations s, these dependencies 
are rather sharp and seem to be mainly determined by the variations of supersaturation 
s(x,T) with x or T. At higher s >0.3, these temperature dependencies become more 
smooth and seem to be determined by an interplay between an enhancement with lowering 
T of both the thermodynamic driving forces promoting the evolution and the vacancy- 
copper atom correlations reducing the copper diffusivity Dq u and thereby slowing down 
the evolution. 

Temporal evolution of the precipitate size distributions Vkif) is illustrated by Fig. 17. 
These distributions are typically rather broad, and they strongly vary with the evolution 
time t. Therefore, the conventional description of these sizes in terms of the mean pre- 
cipitate size R(t) is oversimplified and incomplete. We also find that for the alloy state 
G with a relatively high temperature T=873 K, the size distribution Vf.(t) is much more 
uniform than those observed for the states B and D with the lower temperatures T=773 
K and T=713 K, particularly for the coarsening stage. 

We also describe an improved version of the earlier-suggested stochastic statistical 
method for simulations of precipitation and show that this version can be used for various 
extrapolations of KMC simulations, in particular, for their extensions to the first stages 
of coarsening for which the KMC simulations are time-consuming. 

Comparison of our simulated temporal dependencies for the density and the mean 
size of precipitates in binary Fe-xCu alloys at x=1.34, T=773 K and x=1.15, T=873 
K to the available experimental data [H [21 El E] shows a reasonable agreement within 
both statistical errors of simulations and the scatter of experimental results. The sizes 
N c of critical precipitates calculated by the statistical method of Dobretsov and Vaks [H] 
and presented in Table 1 are close to those estimated in our KMC simulations and in 
experiments by Kolli et al. [5l [6] for the NUCu steels. 

Comparison of results of our simulations for the Fe-1.17Cu and Fe-1.82Cu alloys to the 
data by Kolli et al. [SJ E] about precipitation in NUCu-140 and NUCu-170 steels which 
have the same copper content, xq u = 1.17 and xq u = 1-82, and contain the similar amount 
of other alloying elements enables us to assess the effects of these alloying elements on the 
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precipitation kinetics. The maximum precipitate density d max in each of these two steels 
is lower than that in its binary analogue by about three times. For the nucleation stage, 
the precipitate density d p (t) and their mean size R(t) observed by Kolli et al. j6] in the 
NUCu-140 steel seem to be close to those simulated for the Fe-1.17Cu alloy, contrary to 
the case of the NUCu-170 steel for which the d p (t) and R(t) values at t=0.25 h reported 
by Kolli and Seidman [5] sharply disagree with those simulated for the Fe-1.82Cu alloy. 
In this connection we note that the i?(0.25h)=i?KS value reported in [5] seems to be 
unrealistically large for the early nucleation stage supposed by Kolli and Seidman for the 
NUCu-170 steel at t=0.25 h: i?xs — 4i2 c . Therefore, further experimental studies of the 
nucleation kinetics in Fe-Cu-based steels seem to be desirable. 

For the coarsening stage, the presence of alloying elements in the NUCu steels leads 
to a very strong slowing down of coarsening, by 1-2 orders of magnitude, as compared to 
their binary analogues. A similar strong slowing down of coarsening was also observed 
by Miller et al. jl]) for an irradiated ternary Fe-Cu-Mn alloy. At the same time, for the 
nucleation and growth stages, the effects of alloying elements on the maximum precipitate 
density d max in the NUCu steels and in the Fe-Cu-Mn alloy studied by Shabadi et al. [7J 
seem to be qualitatively different. Some hypotheses about a possible origin of the universal 
slowing down of coarsening in multicomponent Fe-Cu-based alloys with respect to their 
binary analogues are suggested. 
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Research (grant No. 12-02-00093); by the fund for support of leading scientific schools 
of Russia (grant No. NS-215.2012.2); and by the program of Russian university scientific 
potential development (grant No. 2.1.1/4540). 

References 

[1] S.R. Goodman, S.S. Brenner and J.R. Low, Metall. Trans. 4, 2363, 2371 (1973). 

[2] R. Kampmann and R. Wagner, in: Atomic Transport and defects in Metals by Neu- 
tron Scattering, ed. by C. Janotet al., Springer, New York (1986) p. 73. 

[3] M.H. Mathon, A. Barbu, F. Dunstetter, F. Maury, N. Lorenzelli, and C.H. de Novion, 
J. Nucl. Mater. 245, 224 (2004). 

[4] M.K. Miller, B.D. Wirth and G.R. Odette, Mater. Sci. Eng. 353, 133 (2003). 

[5] R.P. Kolli and D.N. Seidman, Acta Mater. 56, 2073 (2008). 

[6] R.P. Kolli, R.M. Wojes, S. Zaucha and D.N. Seidman, Int. J. Mat. Res. 99, 513 
(2008). 

[7] R. Shabadi, R. Taillard, B. Radiget, J. De Baerdemaeker and E. Leunis, Solid State 
Phenomena 172-174, 297 (2011). 

[8] F. Soisson, J. Nucl. Mater. 349, 235 (2006). 

[9] F. Soisson and C.-C. Fu, Phys. Rev. B 76, 214102 (2007). 



17 



[10] F. Soisson and C.-C. Fu, Solid State Phenomena 139, 107 (2008). 

[11] F. Soisson, C Becquart, N. Castin, C. Domain, L. Malerba, and E. Vincent, J. Nucl. 
Mater. 406, 55 (2010). 

[12] A.Yu. Stroev, I.R. Pankratov and V.G. Vaks, Phys. Rev. B 77, 134203 (2008). 

[13] K.Yu. Khromov, F. Soisson, A.Yu. Stroev and V.G. Vaks, JETP 112, 415 (2011). 

[14] VYu. Dobretsov and V.G. Vaks, J. Phys.: Condensed Matter 10, 2261, 2275 (1998). 

[15] V.G. Vaks, I.A. Zhuravlev and K.Yu. Khromov, Zh. Exp. Teor. Fiz. 138, 902 (2010) 
[JETP 111, 796 (2010)]. 

[16] V.G. Vaks and K.Yu. Khromov, Zh. Exp. Teor. Fiz. 133, 115 (2008) [JETP 106, 94 
(2008)]. 



[17] K.Yu. Khromov, V.G. Vaks and I.A. Zhuravlev, http://arxiv.org/abs/ 



V.G. Vaks, LA. Zhuravlev , Zh. Exp. Teor. Fiz. 142 (2012), in press; 
http://arxiv.org/abs/1201.6056l 



[19] M. Nastar, VYu. Dobretsov and G. Martin, Phil. Mag. A 80, 155 (2000). 
[20] F. Soisson and G. Martin, Phys. Rev. B 62, 203 (2000). 

[21] A.V. Barashev, S.I. Golubov, D.J. Bacon, P.E.J. Flewitt, T.A. Lewis, Acta Mater. 
52, 877 (2004). 

[22] I.M. Lifshits, V.V. Slyozov, J. Phys. Chem. Solids 19, 35 (1961). 
[23] C.Z. Wagner, Electrochem. 65, 581 (1961). 



18 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

R (nm) 



10 



26 



10 



25 



10 



24 



10 



23 



I I I 

- 


1 1 

- 

c 






y 

y 






: A 




0'f 


▲ 


1 

1 1 1 









0.5 



1.5 



2.5 



100 



10 



0.1 



I 




i i 


- 




■ 


- 


■ 






Q \ 


• 


I 


A\ 

B ^ 
• 


X A 

c 



0.5 1 1.5 2 2.5 3 

*Cu (%) 



10 



26 



10 



25 



CO 



10 



24 



10 



23 




600 650 



700 



750 

r(K) 



800 850 



900 



C/3 

O 



100 - 



10 - 



1 - 



0.1 - 



0.01 



I 

/Of 


i 


i 


■ 1 


i 


















A . 

/ 

/ 

/ 

/ 

/ _ 

S 

• 

/ 






■ 


4 

/ 


/ 

t 








□ /' 

✓ 

/ 

/ 

A/ 






— -B 








^„A 


4 


--" G : 


A— - 


—A"'" 














A 




G>— 


— e~- - 


— -0 






E 


D 


B 






i 


i 


• 

i 


i 


1 



600 650 700 750 

r(K) 



800 850 900 




50 100 150 



50 100 150 200 250 300 



100 200 300 400 



500 1000 1500 



Bl 



B2 



B3 



B4 



Dl 



D2 



D3 



D4 



20 40 60 80 



25 50 75 100 125 



50 100 150 



100 200 300 400 500 k 




10 



20 



30 



40 



25 50 75 100 125 



50 100 150 200 250 




0.0 1 — — — 1 

10" 2 10" 1 10° 10 1 

t (hours) 




0.0 1 — — — 1 

10" 2 10" 1 10° 10 1 

t (hours) 



